International Journal of Modern Physics A 
© World Scientific Publishing Company 



Anguleir Correlation Estimates for Ultrahigh Energy Cosmic Rays 

CHAD FINLEY* 
for the High Resolution Fly's Eye (HiRes) Collaborationt 

Physics Department, Columbia University, 
538 West 120th Street, New York, NY 10027, USA 



Anisotropy in arrival directions of ultrahigh energy cosmic rays offers the most direct 
way to search for the sources of these particles. We present estimates of the angular 
correlation in the HiRes sample of stereo events above 10^® eV, and in the combined 
sample of HiRes and AGASA events above 4 X 10^^ eV. 
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1. Introduction 

Identifying the sources of ultrahigh energy cosmic rays remains one of the major 
challenges in astrophysics. If these cosmic rays are not too severely deflected by 
magnetic fields while in transit, then their arrival directions will provide crucial 
information about their origins. 

The High Resolution Fly's Eye (HiRes) experiment, consisting of two air fluores- 
cence detector sites, has been making stereo observations of ultrahigh energy cosmic 
rays since 1999. The angular resolution (~ 0.6°) achieved by this experiment sets 
the stage for the most precise measurements yet of small-angle correlation among 
the highest energy cosmic rays. In this paper we analyze HiRes events with energies 
above 10^^ eV observed between 1999 December and 2004 January. More details 
about this data set are given in Ref. ^ We also analyze a combination of events 
from HiRes and from the Akeno Giant Air Shower Array (AGASA) with energies 
above 4 x 10^^ eV. 

2. Analysis 

A standard tool for studying anisotropy is the angular two-point correlation 
function. We define the estimator for the correlation function as w{0) = 
N(9)/ {Nmc{d)) — 1, where N(9) is the number of pairs of events in the data sample 
with angular separation less than 6, and (Nmcid)) is the average number of such 
pairs in simulated isotropic sets with the same number of events and same detector 
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acceptance in right ascension and declination as the data sample. Note that this 
definition of w{9) is cumulative over angles up to 9. It reveals how the correlation 
signal varies as a function of the angular threshold for defining a pair. 

Fig. 1 shows the angular correlation for the HiRes stereo events with energies 
above 10^^ eV, along with the (Poisson) errors y/N / {N,nc} (in general, the variance 
of correlation estimates may be larger than the Poisson variance). No significant 
excess of pairs is observed for any angular scale. 

In 1999 and 2000 (Refs. E| E| the AGASA coUaboration reported significant 
clustering among events above 4 x 10"'^^ eV. In a set of 57 events, seven pairs were 
observed with angular separation less than 2.5°, against an expected {Nmc) = 1-5 
pairs, yielding w(2.5°) — 3.7. The probability for seven or more pairs evaluated 
using Monte Carlo is ~ 0.1%. However, as has been noted in Ref.|3 this probability 
is not indicative of the statistical significance, because of bias introduced by the 
a posteriori choice of the energy threshold and angular separation. 

The angular correlation estimate can be improved by combining the AGASA 
and HiRes data. By itself, the set of 27 HiRes events above 4 x 10^^ eV has no pairs 
separated by less than 5°, yielding w — but with large uncertainty. However, 
the combination of the two sets substantially increases the statistical power over 
either one individually, because of the cross-correlation power of the two sets. 

To determine the expected number of pairs (Nmc), we simulate combined sets 
with 57 events generated using the AGASA acceptance and 27 generated using the 
HiRes acceptance. Although these acceptances are not identical, they have a large 
overlap; the resulting value of (iV^c) is in fact roughly comparable whether one uses 
all AGASA events, all HiRes events, or a combination. 

The two detectors also have different angular resolutions. While this does not 
affect simulated isotropic sets, it could affect the correlation estimate for a real 
clustering signal: a clustering signal of higher significance could appear at smaller 
angular scales because the HiRes angular resolution is several times sharper than 
that of AGASA. However, since 'w{9) includes all pairs with separations less than 
6, this effect will not lead to a reduced signal at larger angles. Therefore rather 
than attempt to estimate the optimal scale for a clustering signal, for the sake of 
comparison we will continue to evaluate ?x;(2.5°) for the combined data set. 

Figure 2 shows the results for the 57 AGASA events alone, and for the 84 
AGASA and HiRes events evaluated jointly. The addition of the HiRes data brings 
one new pair with an AGASA event within 2.5°, yielding = 8, Nmc = 3.0, and 
w(2.5°) = 1.7. The fraction of simulated sets with eight or more pairs is ~ 1%, but 
it must be emphaiszed that this does not represent a chance probability because it 
includes the same bias in the AGASA data set noted above.'' On the other hand, 
an observation of w{2.5°) = 3.7, as seen originally in the AGASA data alone, would 

similar but unbiased analysis using only the AGASA events observed after the choice of cuts 
(i.e. 27 AGASA and 27 HiRes events), results in = 2, N^c = 1.2 and w{2.5°) = 0.7, for a 
chance probability of 34%. Details will be presented in a separate paper. 
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have meant the observation of 14 pairs in the combined data set, corresponding to 
a 10~^ deviation from isotropy. 

3. Conclusions 

Angular correlation estimates for the HiRes events above 10^^ eV reveal no signif- 
icant deviations from isotropy at any small angle. Combining HiRes and AGASA 
events above 4 x lO""^^ eV improves the statistics above this threshold relative to the 
AGASA data set alone, and results in substantially reduced angular correlation. 
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Fig. 1. the HiRes data set above 10^^ eV. Top: Number of observed (A'^) and expected {{Nmc)) 
pairs of events as a function of maximum separation angle 9. Middle: angular correlation w and 
associated uncertainty (see text) for an isotropic set (solid curve) and for the data set (vertical 
bars). Bottom: the fraction of simulated Monte Carlo sets with Nmc > N. Note: values of w in 
neighboring bins are correlated. 
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Fig. 2. the AGASA and the combined HiRes/AGASA data sets above 4 X 10^^ eV. Top: Number 
of observed (N) and expected {{Nmc)) pairs of events as a function of maximum separation angle 
8. Middle: angular correlation w and associated uncertainty (see text) for an isotropic set (solid 
curve) and for the data set (vertical bars). Bottom: the fraction of simulated Monte Carlo sets 
with N„ic > Note: values of w in neighboring bins are correlated. 



